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Goals: 

1.  Identify  all  TOIs. 

2.  Assess  technology. 

3.  Keep  at  least  75% 
non-TIOs  in  ground 


Camp  Beale  TOI 

Main  TOI 


Fuze  and  Fuze-parts  as  TOI 
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Challenges  at  CBE 


□  Magnetic  soil  noise. 

□  Small  S/N  ratio. 

□  TOI  types:  There  were  six  types  of  main  TOIs  and  eight  types  of  native  TOIs. 

□  There  were  only  a  few  (1  to  5)  fuzes  of  the  same  size; 

□  Small  size  (3cm  and  5  cm)  unexpected  native  fuzes  as  TOI. 

□  Multi  Targets 
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UXO  classification 


The  entire  UXO  classification  process  can  be  divided  into  three  parts: 


1 .  Data  Acquisition  2.  Feature  extraction 


P=F-'  [d  ] 
Inverse  Operator 


3.  Decision 


Feature  selection 
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Clustering  and  Classification 
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Forward  Models 


Multi  dipole  mode 


Tx  Rx 


The  scattered  EMI  field  is  approximated 
as  superposition  of  magnetic  fields  from 
each  individual  dipole,  using  the  dyadic 
Green’s  function: 

i= 1 

where 

cf,(r)=  1  ThR,R,  -  7)  ;  R,  =  r,.-r 

4 7tR3  V  / 


The  scattered  EMI  field  is  approximated 
as  magnetic  field  from  groups  of 
interacting  dipoles  using  an  ortho- 
normalized  function  expansion: 

W(v)  =  YJPq(r)-bi, 

q= i 

_  _  q- 1  _  _ 

where  Yq(r)  =  G?(r)  -  £  ¥k(v)  •  Aqk\ 
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Forward  Models 


Multi  dipole  mode 


ONMVS 


Tx 


Rx 


>  nrij  are  determine  from  the  measured 
data  by  solving  a  linear  system  of 
equations. 

>  Uses  individual  dipole 
polarizabilities  for  classification 


>  First  it  determines  bq  from  the  measured 
data  without  solving  a  linear  system  of 
equations,  then  it  backs  out  nij 

>  Uses  total  ONVMS/effective 
polarizabilities  for  classification 
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Inverted  total  ON  VMS  for  different  objects  at 

CBE-MM  data 


Case-138 

io4  . 


Time  [mSec] 


Time  [mSec] 


Time  [mSec] 
Case-464 

104  - T T T - T . . . 


Time  [mSec] 


Case-1572 

104  - T T T - T . . . 


Time  [mSec] 


Case-1534 
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QC  using  Joint  Diagonalization 

Build  a  Multi  Static  Response  (MSR)  data  matrix  [W]  as: 

[ W] (tk )  =  [ht  (tk )]  •  [H(tk )] ,  k=l ,  2, ....  Nt 


Time  Channels  1  :Nt  Represent  MSR  data  matrix  [W]  in  Space  -Time 


[W]  =VD(t|)VT 

•  V  is  a  Eigen 
vectors 
matrix 

•  D(t-|)  is  a 
diagonal 
matrix 


Max(NTx  NRx) 


Max(NTx  NRx) 


The  goal  is  to: 

>  determine  the  eigenvalues  of  [W]  matrix  for  each  time  channel. 

>  find  an  eigenvector  V  that  will  be  shared  by  all  matrices. 


D(tk )  =  VT[W(tt  )]V ,  k=l,  2, 


•  •  • 


t 
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TONVMS  [Arb] 


Training  data  selection 


Choose  discrimination  features 
size  and  decay 


Cluster  the  discrimination  features 


10 

Time  [mSec] 


Identify  an  anomaly  from  each  cluster, 
build  a  custom  training  list. 


r 


W.:0  o 

•  <sr.  ® 

®#  & 


2  3 

Size 
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Camp  Beale  MPV  Classification  Approach 

>  Invert  all  data  for  one,  two  and  three  targets; 

>  Create  custom  training  data  using  a  statistical  clustering  (for  example  Matlab 
function  “clusterdata”,  with  “Ward  Linkage”)  and  time  decay  curves  of  the  inverted 
ONVMS. 

>  Request  Custom  training  data  and  update  TOI  library. 

>  Invert  as  N=4, 5  targets,  if  necessary. 

Three  targets  inversion  for  cell  #758 


<Lio° 

c n 

> 

z 

O 


10 

Time  [mSec] 


Library  target 
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Camp  Beale  MPV  Classification  Approach 


At  CBE  MPV2  suffered  from  magnetic  soil. 
This  was  mostly  due  to  the  deployment 
strategy,  i.e.  the  sensor  was  placed  on  the 
ground. 


Single  object  inversion  results 


The  ONVMS  technique  was  able  to  model  the 
magnetic  soil’s  EMI  responses  as  a  response  from 
a  target.  The  technique  extracted  the  real  target 
features  from  the  noisy  data. 


Multi  objects  inversion  results 


Camp  Beale  MPV  Classification  results 


_fuzesAsClutter 

BealeTrees  Dartmouth  Advanced  Models  None  MPV  Custom  si  vl  UXO 
100 


Fuzes  as  Clutter 


100  200  300  400  500  600  700  800 


Number  of  Clutter  Items  Incorrectly  Classified 
_fuzesAsUXO 


BealeTrees  Dartmouth  AdvancedModels  None  MPV  Custom  s4  v2  UXO 


>  95  anomalies  were  requested  for 
custom  training,  out  of  those  16 
were  TOI  and  79  were  scrap. 

>  All  available  MPV2  data  were 
inverted  and  analyzed. 

>  No  False  Negatives:  all  TOI,  total 
124=  89  (UXO)  +  35(Fuzes),  were 
indentified  correctly. 

>  200  holes  with  clutter  dug, 

>  587  holes  with  clutter  were  not  dug. 

i.e  ~75  %  of  non-TOI  left  in  the 
ground. 
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Eigen  Values  [Arb] 


Camp  Beale  2x2  Array  Classification  Approach 

>  Conduct  data  QC  using  the  JD  technique,  determine  the  number  (N)  of 
potential  targets. 

>lnvert  all  data  as  N  targets,  with  N  is  estimated  from  the  JD. 

>  Create  custom  training  data  using  a  statistical  clustering,  and  eigenvalues  time 
decay  curves. 

>  Request  Custom  training  data  and  update  TOI  library. 


>lnvert  all  data  as  more  than  N  targets,  if  necessary. 


The  JD  eigenvalue  analysis  indicated  that 
cell  #758,  contained  more  than  one 
targets  and  signal  from  a  target  was 
highly  disturbed,  therefore  the  cell  was 
included  in  the  training  data. 


10° 

Time  [sec] 
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Percent  of  UXOs  Correctly  Classified  (%)  Percent  of  UXOs  Correctly  Classified  (%) 


Camp  Beale  2x2  Array  Classification  results 


_fuzesAsClutter 

BealeTrees  Dartmouth  AdvanceModels  None  2x2  Custom  s3  v3  UXO 


_fuzesAsUXO 

BealeTrees  Dartmouth  AdvanceModels  None  2x2  Custom  s5  v3  UXO 


>  98  anomalies  were  requested  for 
custom  training,  out  of  those  24  were 
TOI  and  74  were  scrap. 

>  All  available  2x2  array  data  were 
inverted  and  analyzed. 

>  No  False  Negatives:  all  TOI,  total  of 
124  =  89  (UXO)  +  35(Fuzes),  were 
indentified  correctly. 

>191  holes  with  clutter  dug, 

>  596  holes  with  clutter  were  not  dug. 

i.e  ~76  %  of  non-TOI  left  in  the  ground. 
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Camp  Beale  CH2MHILL  MM  Classification  results 


_fuzesAsClutter 

BealeOpen  Dartmouth  AdvancedModels  None  MetalMapperC  Custom  s3  vl  U 
100 


200  400  600  800  1000  1200 

Number  of  Clutter  Items  Incorrectly  Classified 


fuzesAsUXO 


BealeOpen  Dartmouth  AdvancedModels  None  MetalMapperC  Custom  s3  vl  I 


>132  anomalies  were  requested  for 
custom  training,  out  of  that  25  were  TOI 
and  107  were  scrap. 

>  All  available  CBE  CH2MHILL  data  were 
inverted  and  analyzed. 

>  No  False  Negatives:  all  TOI,  total  170= 
137  (UXO)  +  33(Fuzes),  were 
indentified  correctly. 

>183  Holes  with  clutter  dug, 

>1117  holes  with  clutter  were  not  dug. 
i.e  ~86  %  of  non-TOI  left  in  the  ground. 
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Camp  Beale  Parsons  MM  Classification  results 

Obtained  by  Sky  Research  Production  team: 

Step  1 : 


Case-01616 


Case-01616 


Case-01616 


Case-01616 


105 


Time  [mSec] 
ISO 


Time  [mSec] 


Time  [mSec] 
flare 


Time  [mSec] 


Time  [mSec] 


Time  [mSec] 
37 


Case-01616 


Time  [mSec] 


Time  [mSec] 


■  Used  test  pit  data  collected  at  site  to  establish  feature  libraries  (grey 
lines)  for  potential  TOI 

■  Plotted  anomaly  total  ONVMS  (green,  blue,  black  lines)  against  feature 
libraries  to  identify  high  probability  TOIs  (37mm  shown  in  red  outline) 
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Step  2:  Identify  the  anomalies 


Time  [mSec] 


Fuze 


105 


Time  [mSec] 


Type  2 


Survey  Nail 


Time  [mSec] 


■  Identified  “suspicious”  anomalies  that  did  not  match  feature  library  responses 

■  Requested  ground  truth  for  select  group  of  possible  TOIs 

■  Discovered  two  TOI  fuze  types  and  two  clutter  types  (horseshoe  and  survey  nail) 

■  Added  the  new  TOI  and  clutter  polarizabilities  to  the  feature  library 
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Camp  Beale  Parsons  MM  Classification  results 

obtained  by  Sky  Research  Production  team 


_fuzesAsClutter 

BealeOpen  Dartmouth  AdvancedModels  None  MetalMapperP  Custom  s3  v2  L 


200  400  600  800  1000  1200 

Number  of  Clutter  Items  Incorrectly  Classified 


> 


> 


> 


fuzesAsUXO 


BealeOpen  Dartmouth  AdvancedModels  None  MetalMapperP  Custom  s3  v2  U 


> 

> 


69  anomalies  were  requested  for  custom 
training,  out  of  that  19  were  TOI  and  50 
were  scraps. 

All  available  CBE  Parsons  data  were 
inverted  and  analyzed. 

No  False  Negatives:  all  TOI,  total  170=  137 
(UXO)  +  33(Fuzes),  were  indentified 
correctly. 

253  Holes  with  clutter  dug, 

1047  Holes  with  clutter  were  not  dug. 
i.e  ~81  %  of  non-TOI  left  in  the  ground. 


Targets  that  produced  Unnecessary  digs 


Time  [mSec] 
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Eigen  Values  [Arb] 


Difficult  cases:  JD  applied  to  CBE  anomaly  #2277 


10 

Time  [sec] 


The  curves  of  the  MM  MSR 
data  matrix  eigenvalues  versus 
time  do  not  show  any  evidence 
of  the  37  mm  target. 


4.3312 

4.3312 

CD 

C 

iZ 

t: 

£3312 

4.3312 


The  ground  truth  revealed, 
in  fact,  that  the  37mm 
target  was  90  cm  away 
from  the  MM  center. 


xIO 


MM  Position 


37  mm 


Barbered  wires 

•  •  % 


EM61  flag 


6.4711  6.4711  6.4711  6.4711  6.4711 


Easting 


x  ICf 
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Summary 

>  Advanced  EMI  models  were  applied  to  CBE  Cued 
Data  sets. 

>  Studies  showed  that  the  models  are  able  to  deal 
effectively  with  cluttered  environment. 

>  Classifications  were  demonstrated  for  both  R&D  and 
production  EMI  sensors. 

>  No  False  Alarms. 

>  The  models  are  easy  to  use  for  general  users. 

>  The  models  were  able  to  classify  targets  as  small  as 
3  cm  size  fuzes. 
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